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SUMMARY 


A numerical lifting-surface method has heen used to calculate direct gust 
forces and moments on wings of several planforms. The gust velocities are con- 
tinuous and vary sinusoidally in the stream direction and are also uniform across 
the wing span. The procedure has the advantage of rapid machine calculation and 
includes the effects of wing planform, nonsteady subsonic flow, and induced flow 
effects. The method provides for calculation of gust forces on a basis consis- 
tent with that for the calculation of forces due to motion and deformation. The 
results include the in-phase and quadrature components of the following 
quantities: (a) spanwise distribution of section lift coefficient, (b) total 

lift coefficient, and (c) total pitching-moment coefficient. In addition, gen- 
eralized gust forces on approximate fundamental cantilever bending modes (para- 
bolic) are also included. Results have been obtained for 60° and 75° delta wings, 
55 ° sweptback wings of aspect ratio 4.00 and a 5° sweptback wing of aspect 

ratio 11.60, and an unswept wing of aspect ratio 6.00. Conditions for which cal- 
culations were made include two Mach numbers (generally, 0.40 and O. 9 O) and a 
reduced-frequency range of 0 to 1.0. The direct gust forces and moments are in 
forms suitable to be inserted in equations of motion used in the calculation of 
the dynamic responses of flexible lifting vehicles to random turbulence and to be 
compared with results from other methods. 


INTRODUCTION 


Current analytic methods for calculating airplane responses to continuous 
vertical atmospheric turbulence (random process theory) require a knowledge of 
unsteady generalized aerodynamic forces due directly to turbulence velocities. 
Within the framework of random process theory, these aerodynamic forces are 
described in terms of the forces produced directly by vertical gust velocities 



which vary sinusoidally along the flight path of the airplane and which are uni- 
form along the span of the lifting surface. For convenience in application, 
values of unsteady aerodynamic forces due to these inputs are desired for the 
planform of interest. 

Explicit solutions of the aerodynamic equations for a general lifting surface 
have not been found. Approximations involving strip theory have been utilized, 
such as the use of steady- state spanwise lift distributions together with an 
amplitude function based on two-dimensional flow and simple time lags to account 
for gust field penetration. Results of work on unsteady aerodynamic forces perti- 
nent to gust analysis can be found in references 1 to k. Information for wings 
of finite aspect ratio in subsonic compressible flow is very limited. 

Reference 5 has provided a numerical method for solving the equations for a 
lifting surface (utilizing a high-speed digital computer) which can be used for 
calculating sinusoidal gust forces in subsonic flow. This method was developed 
primarily for flutter studies, but the applied downwash description permits it to 
be extended to the gust case. It therefore provides aerodynamic forces on a con- 
sistent basis for both oscillating airstream and oscillating lifting surfaces. 

An application of this technique to a swept wing has been demonstrated as an inci- 
dental part of the study of reference 6 with satisfactory results indicated. 

The purpose of the present report is (l) to demonstrate the application of 
the procedure of reference 5 to the calculation of generalized forces, including 
lift and pitching moments, due to sinusoidal gusts for several commonly used wing 
planforms and (2) to present these results in the form of curves that may be uti- 
lized in simplified calculations. The data for planforms considered herein 
include results at two subsonic Mach numbers for 60° and 75° delta wings, 35° 
sweptback wings of two aspect ratios, one 5° sweptback wing, and two unswept wings. 
A comparison is made of the results of the numerical procedure and those of a 
modified strip theory. These results will add to the limited available Informa- 
tion on unsteady lift functions for some typical wings of finite aspect ratio in 
sinusoidal gusts in subsonic flow. 


SYMBOLS 


A 

Af,Bi 


* * 

Ai,Bi 


b 

C i> D i 


aspect ratio, b^/s 

real and imaginary parts, respectively, of spanwise distribution of 
generalized gust force per unit gust velocity 

normalized values of A^ and B^, respectively (see section entitled 
"Gust Forces") 

wing span, ft 

real and imaginary parts, respectively, of total generalized gust force 
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C* p* 


C L 

% 


c a 

c l 

c r 

c t 

K 

k 

L 

M 

Mij 


My 


p,Q 


Ap 

*5 

Ap G 


Qi(“) 

Q? 


normalized values of and T>±, respectively (see section entitled 

"Gust Forces") 

lift coefficient, L/qS 
pitching-moment coefficient, 2My/qSc r 
lift-curve slope, dCjj/da, per radian 

local wing chord, ft 
average chord, S/b, ft 

section lift coefficient 

root chord of lifting surface, ft 

wing chord at tip, ft 

kernel of equation (2), integral equation (see ref. 5) 
reduced- frequency parameter, mc r /2U 
lift-force, lb 
Mach number 

generalized inertia force or moment coefficient representing coupling 
between ith and jth modes, slugs 

pitching moment about midpoint of root chord, positive nose up 

real and imaginary parts, respectively, of normalized unsteady lift func- 
tion for wing of infinite aspect ratio; for steady flow: P = 1, Q = 0 

differential pressure, lb/sq ft 

pressure distribution due to downwash of jth mode, per unit generalized 
coordinate, ^Ap^Sqj 

pressure distribution due to gust downwash, per unit downwash velocity 

ith mode generalized aerodynamic force or moment 
generalized aerodynamic force due to unit gust velocity 
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q1 (oj) generalized aerodynamic force due to wing motion 
q dynamic pressure, pU^/2, lb/sq ft 

q. ((Jo),q i (u>) ith and jth generalized coordinates, respectively 


s 

t 

u 

w 

x 

y 


wing area, sq ft 
time, sec 

velocity of airstream, ft/sec 
applied dovnwash velocity, ft/sec 
gust dovnwash velocity, ft/sec 

coordinate along chord of lifting surface measured rearward from 
leading edge of root chord 

coordinate along semi span of lifting surface measured outward from 
root chord 


y 


* = 


y 

b/2 


z wing deflection, ft 

a angle of attack, radians 

A sweepback angle of l/4 chord line of lifting surface, deg 

Ai e sweepback angle of leading edge of lifting surface, deg 


A 

4i(x,y) 


wavelength, ft /cycle 
taper ratio, c-t/ c r 

shape of ith mode of airplane structure 


p density of air, slugs/cu ft 

a) frequency, radians/sec 

Subscripts: 

i,j integers indicating modes 

A bar above a symbol indicates amplitude only. 
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ANALYSIS 


In this section the relationship of the gust force terms to the other terms 
in the general equations of motion will he shown. The numerical procedure for 
determining the pressure distribution over the lifting surface will also he 
briefly described. 


Gust Forces 

The general equations of motion based on uncoupled modes of motion and defor- 
mation for an airplane flying through a continuous sinusoidal gust field can be 
written as follows (based on refs. 6 and 7 ): 


n 


-^ 2 M ij<lj(^) + o^?M ii (cD)q i ((n) = Q i (aj) « YJ % jfatejfa) + QjfaW 

1 J=1 


(i — 1, 2 >3,? • • • n ) (l&) 

where the terms Q^j(to) and Q^(ao) are the generalized aerodynamic forces due 

to airplane motion and gust velocity, respectively. The generalized forces on 
the right-hand side of the equation are defined as 


and 


= JJ !i(x,y) (x,y;u>) dx dy 

S j 


Q?(<o) = JJ li(x,y) (x,y;o>) dx dy 


(lb) 


(lc) 


This report is concerned with the determination of generalized gust forces 
Q!?(o>) . Integrating equation (lc) over the chord yields 


Q?M 



(Id) 


where and represent the real and imaginary parts, respectively, of the 

spanwise distribution of the generalized gust forces per unit amplitude of gust 
velocity. The terminology used herein of real and imaginary implies the in-phase 
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and out-of -phase components of a quantity with respect to the applied downwash 
at the leading edge of the root chord of the lifting surface. Integrating equa- 
tion (id) over the span yields 


Q^(cd) = Cp(cu) + iD-j^co) (le) 

where Cp and Dp represent the real and imaginary parts, respectively, of the 
total generalized gust forces. 

For man y applications the degrees of freedom of rigid-body translation, 
rigid-body pitch, and first flexible bending mode of vibration are sufficient to 
determine the desired result. Consequently, the results presented herein are 
responses for these degrees of freedom. These results are as follows, referring 
to equation (lc): 

(1) For rigid-body translation, i = 1 and |p = lj therefore, Ap and Bp 
may be interpreted as the real and imaginary parts of the spanwise lift distri- 
bution. The terms Cp and Dp are then the real and imaginary parts of the 
total lift obtained by integrating Ap and Bp over the span of the lifting 
surface. For purposes of emphasizing the nonsteady flow characteristics, norma- 
lized values of Cp and Dp are obtained from in tegration after dividing Ap 

and Bp by the modulus of the steady-state lift ifCp + Dp. These normalized 

* * V 

values will be designated Cp and Bp. 

( 2 ) For purposes of presentation, the quantities Ap and Bp have been 
normalized by dividing each value of Ap and Bp by the magnitude of the total 

lift for the frequency under consideration. This procedure eliminates the 
effect of magnitude and illustrates the effect of reduced frequency on the 
in-phase and quadrature portions of the spanwise lift distribution. These nor- 
malized values will be designated A* and Bp. 

(3) For rigid-body pitch, i = 2 and thus, C2 and B2 may be 

interpreted as the real and imaginary parts of the total pitching moment. As 
was the case for Cp and Bp, C2 and D2 will be presented as normalized 

values C* and These are obtained by dividing A2 and B2 by the modulus 

r 2 2 

of the steady- state moment 1/ + I>2* 

( 4 ) For a first bending mode approximated by a parabolic spanwise shape ; 

i = 3 and | ^ = (y*) 2 . Since | ^ is only a function of y, A^ = Ap(y *) 2 and 

Bj = Bp(y*) 2 . Spanwise integration yields Cj and which are the real and 

imaginary parts of the gust forces for this mode shape. 
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In this paper A* and B*, C* and D*, C* and I>2 > and Cj/qS and Dj^qS 
will subsequently be presented for the various planforms considered. 


Determination of Pressure Distribution 

The real and imaginary components of the pressure distribution Ap G / w g in 

equation (lc) are obtained by the numerical method of reference 5 with the sub- 
stitution of an expression for the gust downwash for that of motion and deforma- 
tion dovnvash. A brief description of some of the features of the numerical 
method and of the expression for the gust downwash is given below. 

The numerical method is applied to the solution of a linear integral equa- 
tion of the form 

w(x,y) = JJ' | , t])k(x-| , y-r|jk,M) d| dT) (2) 

q S 

The function K, the kernel of the integral equation, physically represents the 
contribution to the downwash at a field point (x,y) due to a sinusoidally pul- 
sating pressure doublet of unit strength located at any point The solu- 

tion of the equation is obtained from a high-speed digital computer program which 
requires only that the following information be provided: 

(1) The coordinates of nine- 1 - control points on one-half of the wing and the 

downwash at each point 

(2) The parameters that describe the wing boundaries 

(3) The subsonic Mach number 

(4) The values of the reduced frequency k = o>c r /2U for which solutions 

are desired 

The method of reference 5 was developed for flutter studies; consequently, 
the expression given therein for the downwash is in terms of sinusoidal wing 
motion and deformation. This expression can be written 


w(x,y)e :kr>t = 


u|£(x,y) + iuz(x,y) 
ox 


e tot 


(5) 


For the generalized gust forces no motion or deformation of the wing is 
desired in the expression for the downwash. Rather, the downwash expression must 
account for the sinusoidal variation of the gust downwash along the wing in the 
chordwise direction. The appropriate expression is 


ISince preparation of this report, the computer program has been extended to 
12 and l6 control points. 
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w(x,y)e icDt 


(4) 



where x c is the location of a control station relative to the leading edge of 
the root chord. Equation (4) can he written 

w(x,y)e iart = Wge iayt fcos 2k— - i sin 2k— J (5) 

& \ c r c r / 


where 


k 


03C r 

“2U 


For this downwash the pressure distributions are obtained numerically and the 
generalized gust force is obtained by integrating over the lifting surface. (See 
ref. 5*) 


RESULTS AMD DISCUSSION 


Presentation of Results 

The planforms for which results are presented herein are shown in figure 1. 
These planforms include 60° and 75° delta wings and two wings with quarter- chord 
sweepback angle of 55° and taper ratio of 0.42, with aspect ratios of 9-^3 and 
4.00. Also included is a wing with a 5° sweepback angle, aspect ratio of 11.60, 
and taper ratio of 0.44. Unswept wings are also considered; the results for the 
"classic" rectangular wing of aspect ratio 6,00 are included for comparison 
along with an unswept wing (quarter -chord line perpendicular to root chord) of 
aspect ratio 6.00 and taper ratio 0 . 50 . 

Calculated values of steady-state lift- and pitching-moment -curve slopes for 
all planforms considered are shown in table I. Some of these values are discussed 
in a subsequent section entitled "Limitations." 

Values of the normalized real and imaginary parts of the generalized gust 
force for vertical translation (total lift), Cq and Dq, respectively, are 
given in figures 2(a) to 2(h) . The phase relationship is given with respect to 
the downwash at the leading edge of the root chord. Figure 2(a) gives results 
from reference 1 (generally termed P and Q) which are for a wing with infinite 
aspect ratio in incompressible flow (M - 0) and provides a convenient reference 
for the results for all the other planforms. Figure 2(b) gives results for a 
rectangular wing of aspect ratio 6.00 in incompressible flow. A comparison of 
present results with those obtained by using the approximation given in equa- 
tion (33) of reference 4 is shown in figure 2(b) and shows almost the same 
results . 
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Values of the normalized real and imaginary parts of the generalized gust 
force for rigid pitching motion (total pitching moment) are given in figures 3(a) 
to 3(f) as C2 and Dp- This pitching moment is about the midpoint of the root 

c r 

chord of the lifting surface. The moment relation used is My = C^qS-^-. 

The normalized real and imaginary parts of the spanwise distribution of gen- 
eralized gust forces for rigid-body translation (lift) A* and are shown in 

figures Ha) to 4(m) for the various planforms at two Mach numbers and various 
values of k (generally 0, 0.5, and 1.0). It should be noted again that all 
phase relationships are with respect to the downwash at the leading edge of the 
root chord of each planform. All curves presented would, of course, exhibit sig- 
nificant changes if the downwash reference were shifted to another location. 

Figures 5(a) to 5(g) present values for the real and imaginary parts of the 
generalized gust-force coefficients for a deflection mode with a parabolic shape 
approximating the first wing bending mode Cjj qS and D^/qS. 

Discussion 

Li mi tat ions .- The present procedure provides for the calculation of genera- 
lized gust forces on a basis consistent with that for the calculation of genera- 
lized aerodynamic forces due to motion and deformation of a lifting surface. It 
is the purpose of this section to indicate the range of configurations and condi- 
tions for which the preceding statement applies. The principal limitations are 
with respect to aspect ratio and reduced frequency with a slight limitation on the 
subsonic Mach number range, the maximum Mach number being approximately 0.95- 

The application of the present procedure is limited to the lower values of 
aspect ratio - say to less than 10 or 12. In general, errors increase with an 
increase in aspect ratio. This error as displayed by the steady-state lift-curve 
slope for a rectangular wing in incompressible flow (M = 0) is illustrated in the 
following table. Values calculated by the present procedure are given in the 
second column. Particularly noteworthy is the value for aspect ratio 20 of more 
than 2rt (the limit value for infinite aspect ratio) . 


Aspect 

ratio 

Steady-state lift-curve slope obtained from - 

1 

Present 

method 

Lifting-line 
theory 
(ref. 8) 

Lifting- surface 
theory 
(ref. 9) 

6 

4.28 

4.18 

4.30 

10 

5.34 

a 4.8l 


12 

a 5-68 

5.00 

— 

15 

6.15 


— 

20 

6 .84 


— 


interpolated values.. 


9 




A further indication of the error in the present procedure is provided hy a com- 
parison of the results with those from steady-state lifting-line theory presented 
in the third column. The latter results were obtained from information in refer- 
ence 8 and are considered to be progressively more accurate with increasing 
aspect ratio. In contrast, the results of the present numerical procedure become 
more accurate with decreasing aspect ratio. The tabulation indicates that errors 
in the results of the present procedure are less than about 10 percent for aspect 
ratios less than 10. The source of the errors at the higher aspect ratios is 
thought to be the limited number of control points and the particular numerical 
integration employed. 

An additional assessment of the present procedure may be made from the pre- 
ceding tabulation. The lift-curve slope for an aspect ratio of 6, obtained from 
lifting- surface equations by a different method of solution described in refer- 
ence 9, is given in the fourth column. The close agreement with results from the 
present method indicates consistency of results from applying lifting- surface 
theories at the lower aspect ratios. 

Inasmuch as the aspect ratios of two of the several wing configurations con- 
sidered herein were in the vicinity of 10, the errors with respect to the results 
of lifting-line theory were investigated in a manner similar to that for the 
rectangular wing. The data from lifting-line theory were obtained by interpola- 
tion of results from references 8 and 10. Prandtl-Glauert corrections for com- 
pressible flow were made in accordance with the procedure described in refer- 
ence 10. The results are presented in the following table: 


Planf orm 

A 


M 

Steady-state lift -curve slope 
obtained from - 

Present 

method 

References 8 
and 10 

5° swept 

11.60 

0.44 

0.33 

6 .08 

5.36 

5° swept 

11.60 

.1* 

•TO 

7-55 

6 .65 

35° swept 

9M 

.42 

.58 

5.06 

4.66 

35° swept 

9M 

.42 

1 

•90 

6.32 

5.85 


The variations of error with aspect ratio are consistent with those for the 
rectangular wing. 

There is a value of reduced frequency beyond which errors in the present 
procedure become excessive. This limitation arises because, in the present usage, 
the applied downwash is satisfied at only three points along the chord, at each 
of three semispan stations. In numerical operations it is common practice to 
require at least nine or ten uniformly spaced points to describe a cycle of a 
sine function. Therefore, in the present procedure, errors may be expected to 
increase sharply for frequencies greater than those having a wavelength of about 
3 chords. For all wings considered herein, the maximum chord is the root chord 
and the frequency limitation can be expressed as 
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(6a) 


~h > 3c r 


and 

* = MI = ^ (6b) 

co k-^ 


Therefore, 


^ > 3c r (6c) 

k l 


and 


k-, < | (6a) 

x 3 

or approximately, 

k x < 1 (6e) 


A similar limit in reduced frequency arises for sweptback wings due to a 
sinusoidal variation of downwash along a line through the most forward control 
points (in this case, the quarter-chord line) of the wing and replaces the limit 
described previously for the condition ^^2 < ' ^ an ^ as indicated. The 

following sketch is presented to clarify the terminology. 
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To approximate a sinusoid along the span 


hut 


A = — E > 5 1 sin A 
k 2 


(7a) 


7 


= 'b/g- 

cos A 


= Wing length along quarter-chord line 


(7b) 


hence. 


* 2 < $ (7c) 

0 | tan A 

For the limitation of equation (7c) to apply, 

kg < kp 

and, for this condition, considering equations (6d) and 7 ( c ) > 


rt c r < * 

5 & tan A 5 
2 


therefore , 


-^£- < tan A 

*72 


(7d) 


(7e) 


For example, for A = 35° and 


= 


3«35j it can he seen that 


3-35 


< 0.700j 


or 


therefore, kg is the l imi ting reduced frequency and kg < & — ^ ^ ^00) 

kg < 0 . kk6 . Thus, for this planform the limitation on kg due to sweep is quite 
severe. For numerical procedures using more control points, limitations on kp 
kg may he relaxed. 


or 


Application of tables and figures .- As mentioned previously, some of the 
re suit s were normalized to illustrate better the effects of nonsteady flow in 
terms of reduced frequency. For application to airplane -response calculations, 
these results must he converted to forces, as defined by equation (id), suitable 
for insertion into equation (l) . The necessary conversions are described in the 
following paragraphs. 
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The results in figures 2 and 3 are converted by use of the equations 


Q l = C 1 (k) + iD x (k) = qs[c*(k) + iD*(kj] (8) 

<4 = C 2 (k) + iD 2 (k) = q (9) 


respectively. 

The results in figure 4 are converted to spanwise force distributions suit- 
able for computing generalized gust forces on flexible wing bending modes by use 
of the equation 


A x (y*,k) + lB 1 (y*^k) 



k=0 


J[cJ(k)] 2 + [pjCkj] 2 [Aj(y*,k) + iBj(y*,kj) 

( 10 ) 


The results in figure 5 may Be converted to generalized gust forces on 
parabolic wing bending modes by multiplication by qS. 

The values of and (— £] used in equations (8) to (10) are 

V la, / k=0 \ da, / k=Q 

given in table I. 

Compa r isons with modified strip theory .- In the absence of a procedure or 
the facilities for a rapid numerical solution of lifting-surface equations, there 
is a temptation to consider some form of strip theory. In this section a com- 
parison is made of results of a strip theory with those of the numerical lifting- 
surface calculations . 

One form of strip theory which appears to offer a useful approximation for 
gust-force calculations for wings of moderate to large aspect ratio consists of 
assuming that the lift at a local point on the span is equal to the lift at that 

( c l c 

point based on the steady- state lift distribution 

three-dimensional surface in a uniform downwash (see refs. 8 and 10) multiplied 
by the local angle of attack at the leading edge w i e pj and by the nonsteady 

lift function for a wing of infinite aspect ratio in a sinusoidal gust field 
P + iQ (see ref. 1, for example). 


for the 


'a=Constant 
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This may he expressed as 


where 


C Z C > 

\ w g 

CL c a 

U 


^icjot 


VL\ 

F' a / a=C on s t ant 


w 2e 

U 




w- 


Ze 



tan A Ze\ 
U / 


( 11 ) 


and describes the lag of penetration of the leading edge at the local chord rela- 
tive to the leading edge at the root chord. 

n 

This equation is related to the generalized gust force, for example, by 




Ql = 


<1 SC L, 


'a 


U 


r 

n 


c l c 


0 C L c a 


ay 


or 


^ r 0 fe) 

0 \ L a /<x=Constant 


p (t) + iQ ( k d 


cos j oo— - tan A 


Ze. 


i sin (cjoi tan 

U 


A Ze) 


ay 


(12) 


from which 



where 
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and 


(15) 



- (l - >r)y* 


For purposes of comparison with the results of the numerical lifting- surface 
calculations, A* and Bp from the strip theory were calculated for a 35° 

sweptback wing and for a 60° delta wing. The Prandtl-Glauert compressibility 
correction as described in reference 10 was applied and the nearest geometric 
configuration of sweep, taper ratio, and aspect ratio was used to obtain 

( — \ from reference 8 . The function P(k) + iQ(k) was obtained 

\ 1 "° a / a =c on s t ant 

from results given in reference 1. As indicated in equation (12), in the appli- 
cation the reduced frequency k was based on the local chord. A check has shown 
very little difference between these results and those for P(k) + iQ(k) with 
the reduced frequency based on the root chord. 

The results of the strip theory in the form of Ap and Bp are given, 

together with those of the numerical lifting- surface calculations, in fig- 
ures 6 (a) and 6 (b) for the 35 ° sweptback wing and in figure 6 (c) for the 
60° delta wing. The two sets of results agree fairly well for both wings at the 
frequencies considered and should compare similarly for other frequencies lower 
than the highest shown. Although not presented in the figures, values of C-^ 

and Dp from strip theory are generally somewhat lower than those from the 

lifting- surface calculations. The ratios of amplitudes are about 78 percent for 
the rectangular wing with A = 6.00; 88 percent for the 35° ^wept wing with 
A = 9.^3; and 85 percent for the 60° delta wing. 


CONCLUDING REMAEKS 


The responses of six wing planforms to continuous sinusoidal gusts have been 
calculated by the use of a lifting-surface procedure. This procedure provides 
for the calculation of generalized gust forces on a basis consistent with that 
for the calculation of generalized aerodynamic forces due to motion and deforma- 
tion of a lifting surface. There are known limitations to the procedure, however, 
and these are discussed; the procedure is limited to aspect ratios less than 
about 10. The gust forces and moments are in forms suitable to be inserted in 
equations of motion for calculation of dynamic responses of lifting surfaces to 
random turbulence . 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 20, 1962 . 
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TABLE I.- STEADY-STATE VALUES OE LIFT- AND PITCHING-MOMENT - CURVE SLOPES 


Planf orm 

A 

' Kj , 

M 

dC L /dcx 

dC m/ da 

Unsvept 

6.oo 

1.00 

0 

4.28 

2.49 

Unsvept 

6.00 

.50 

.40 

4.77 

Unsvept 

6.oo 

.50 

.90 

6.90 

3.80 

5° svept 

11.60 

.44 

• 33 

6.08 

1.17 

5° svept 

11.60 

.44 

.70 

7.55 

1.42 

55° svept 

4.oo 

.42 

.40 

3.74 

-1.35 

35° svept 

4.oo 

.42 

.90 

4.78 

-I.78 

35° svept 

9.43 

.42 

.58 

5.06 

-7.94 

35° svept 

9.43 

.42 

.90 

6.32 

-10.17 

60° delta 

2.30 

0 

.40 

2.63 

-.42 

60° delta 

2.30 

0 

.90 

3.07 

-.63 

75° delta 

1.07 

0 

.40 

1.46 

-.30 

75° delta 

1.07 

0 

.90 

1.57 

-.40 
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(a) The two-dimensional wing (A = oo) in incompressible flow (M = 0) . 

and D* are the P and Q functions, data given in ref. l.j 

Figure 2.- The variation with reduced frequency k of the in-phase (real) 
and quadrature (imaginary) components of the generalized gust force 
tending to produce vertical translation (equivalent to total lift), 
normalized by the force magnitude for k = 0, for the various wing 
planforms shown in figure 1 at Mach number M in a continuous sinus- 
oidal gust field. The phase referral is to the gust velocity at the 
wing apex. 
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Reduced frequency, k 


(b) Unswept wing; A = 6.00; Aji = 1.00; M 


Figure 2.- Continued. 





ralized gust force, C- and 





Generalized gust force, C and 



(d) 5° swept wing; A = 11.60; = O.Uk 

Figure 2.- Continued. 
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’orce, C, and 






(f) 35° swept wlngj A = 9-^3} = O.k-2. 


Figure 2.- Continued. 
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Reduced frequency, k 


(g) 60° delta ving; A = 2 ,30; Tsrp = 0 

Figure 2.- Continued. 




Generalized gust force, C-, and 


1.2 


.8 


♦ iH 
« 



o 


-.4 


-.8 



(h) 75° delta wing; A = 1.07; 7^=0. 

Figure 2.- Concluded. 



Generalized gust force 


0.40 

.90 


mi 



Reduced frequency, k 


(a) Unswept wing; A = 6.00; ~hj< = 0*50. 

Figure 5.- The variation with reduced frequency k of the in-phase and 
quadrature components of the generalized gust force tending to produce 
a displacement in pitch about an axis through the root midchord (equiv 
alent to total pitching moment), normalized by the force magnitude for 
k = 0, for the wings of figure 1 in a continuous sinusoidal gust field 
The phase referral is to the gust velocity at the wing apex. 




Generalized gust force, CL and 



(b) 5° swept wing; A = 11 .60; = O.Ml-. 

Figure 3,- Continued. 
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Genera 






Generalized gust force, CL and 


1.2 



(e) 60° delta wing; A = 2.30; Trp = 0. 

Figure 3»- Continued. 
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Generalized gust force, C P and 



(f) 75° delta wing; A = 1.07; Aj = 0. 

Figure 5*~ Concluded. 
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Semispan location, 


(a) Unswept wing; A = 6.00; = 1.00; M = 0. 

Figure 4.- Spanwise distribution of the in-phase and quadrature components of the 
generalized gust force tending to produce vertical translation (equivalent to 
s ection l ift), normalized by division by the magnitude of the total lift 

I / C? + D?, at each value of k for a continuous sinusoidal gust field. 




Generalized gust force, A n and 
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(b) Unswept wing; A = 6.00; Tvy = 0 . 50 ; M = 0.4 
Figure 4.- Continued. 



neralized gust force, A 1 and 




Generalized gust force, A, and 
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(d) 5° swept wing; A = 11. 60; Tvp = 0.44; M = 0.33 

Figure 4.- Continued. 




Generalized gust force } A n and 



(e) 5° swept wing; A = 11. 60; >rji = 0.44; M = O.'JO. 

Figure 4.- Continued. 
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Generalized gust force, A, and 



(f) 55° swept wing; A = 4.00; = 0.4-2; M = 0.40 

Figure 4.- Continued. 
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Generalized gust force, A 1 and 






Generalized gust force, A n and 



k2 


(j) 60° delta wing; A = 2. 30; Typ = 0; M = 0.40. 

Figure h.- Continued. 









Semi span location 




b/2 


(l) 75° delta wing; A = 1.072; 7vp = 0; M 


Figure 4 .- Continued. 




Semi span location, 


(m) 75 ° delta wing; A = 1.072; Arp 


b/2 

= 0; M 


Figure 4 .- Concluded. 














Reduced frequency, k 


(b) Unswept wing; A = 6.00 j >rp = 0.50. 
Figure 5 .- Continued. 




(c) 5° swept wing; A = 11.60} >»ji = 0.44. 

Figure 5 .- Continued. 



Reduced frequency, k 


(d) 35° swept wing; A = 4.00; Tvp = 0.42 

Figure 5*- Continued. 





(e) 55° swept wing; A = 9*^3) = 0.42. 

Figure 5* - Continued. 
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Reduced frequency, k 


(f) 60° delta wing; A = 2.3°; = 0 

Figure 5*- Continued. 





(g) 75° delta wing; A = 1.07; ?*j» = 0. 

Figure 5 .- Concluded. 
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(a) 35° swept wing; A = 9*^3j = 0.42; M = O. 58 ; k = 0.034. 

Figure 6.- A comparison of results from lifting-surface theory and from modified 
strip theory for spanwise distributions of the in-phase and quadrature compo- 
nents of the generalized gust force tending to produce vertical translation 
(equivalent to sectio n lift), normalized by division by the magnitude of the 

total lift Jc? + D? . 




Generalized gust force, Ai and 



(b) 35° swept wing; A = 9*^3; ^ = 0.b2; M = O. 58 ; k = 0.3^3- 

Figure 6.- Continued. 
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Generalized gust force, A-, and 



(c) 60° delta wing; A = 2. 50; Tvp = 0; M = 0.40; k = 0.50. 

Figure 6.- Concluded. 
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